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^ ■ Abstract 

o 

We consider the classical equations of the gravitating Abelian-Higgs 
■ model in an axially symmetric ansatz. More properties of the solutions of 

pH I these equations (the Melvin and the sting branches) are presented. These 

solutions are also constructed for winding numbers N = 2. It is shown 
^ I that these vortices exist in attractive and repulsive phases, separated by 

the value of the Higgs coupling constant parameter leading to self-dual 
equations. 



1 Introduction 



Most of the well known classical solutions in field theory are constructed in 
Minkowski space. Taking into account gravity, one recovers them as limiting 
cases but some new solutions appear which cannot be reached continuously from 
the existing one in flat space. Such a situation occurs when one considers an 
Abelian-Higgs system. Imposing cylindrical symmetry,one finds that there is 
only one solution : the cosmic string . When gravity sets in, this configuration 
experiences a deformation which manifests itself as a conical singularity. As the 
angular deficit vanishes with the coupling constant measuring the gravitation- 
nal interaction, this solution evolves smoothly from fiat to curved space. It has 
been shown recently that there exists another solution exhibiting axial symme- 
try in the gravitating Abelian-Higgs system The profiles of axial functions 
parametrizing the metric and the matter fields have been obtained for a unit 
winding number. 

In this work, we first address the question of the link between the masses 
(inertial mass and Tolman mass) of the new configurations and the parameters 
of the model (the ratio a of the vector mass by the scalar mass on one side, the 
product 7 of the Newton constant by the square of the Higgs-field's expectation 
value on the other side). They are contrasted with those of the corresponding 
gravitating string. We then turn to configurations with higher winding numbers 
and analyze their stability. The phenomenon first pointed out in [0 is present 
for the gravitating strings and for the new (Melvin-like) solution. 

In a certain domain of the parameter space, the solution becomes closed i.e. 
defined within a maximal radius of the axial variable. We also investigate the 
dependence of this maximal radius versus a and 7. 
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2 The equations. 



We consider the gravitating Abelian Higgs model in four dimensions. It is de- 
scribed by the action 

S = f d^Xy/^C (1) 



C = lD,<l>D^r - \f,^F^'^ - - + (2) 

where = V^ — ieA^ is the gauge covariant derivative, is the gauge potential 
of the U(l) gauge symmetry, F^j, is the corresponding field strength and is a 
complex scalar field with vacuum expectation value v. The usual Ricci scalar is 
denoted R. We use the same notations as [jl|. 

In the following we study the classical equations of the above field theory 
within the cylindrically symmetric ansatz. The metric and matter fields are 
parametrized in terms of four functions of the cylindrical radial variable r : 

ds^ = N\r)de - dr^ - L\r)d(t)^ - N\r)dz^ (3) 

= t;/(r)e^"^ (4) 

Af.dx^' = \n - P{r))d(j) (5) 

Here N, L, P, f are the four radial functions of r, n is an integer indexing the 
vorticity of the Higgs field around the axis. In [|l| it was shown that the above 
ansatz consistently leads to the following system of equations for the functions 
N,L,PJ: 



(6) 



P'^ 1 
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where a = e^/A and 7 = StiGv"^ are the physical constants. The dimensionless 
coordinate x = V Af ^ has been used and the primes mean derivatives with respect 
to X. The problem posed by the differential equations is then further specified 
by the following set of boundary conditions 

N{0) = 1 , A^'(O) = (10) 
L(0) = , L'(0) = 1 (11) 
P(0) = n , lim P(x) = (12) 

X— »oo 

/(O) = , ]im/(x) = l (13) 

which are necessary to guarantee the regularity of the configuration at the origin 
and the finiteness of the inertial mass defined next. 

The different solutions of the system can be characterized by the inertial mass 
per unit length 

Min = J ^0 dXidX2 (14) 

and by the Tolman mass per unit length 

Mto = J Tj;^ dxidx2 (15) 

In the cylindrically symmetric ansatz used and with the dimensionless variable 
X = \^Xvr the two quantities above are given by the following radial integrals 

GMin = ^Min 
o 

- if + ^ + ^ + (16) 

GMto = ^Mto 

= 1 rdxN^L(^--{l-f)A (17) 

The integral determining the Tolman mass can be evaluated by means of 
equation (^ and is given by the following limit 

GMto = l lim {LNN') (18) 
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3 Discussion of the global solutions. 

Let us first discuss the solutions of the system @, (P) which are regular on [0, oo], 
i.e. when N, L have no zero on this interval. 

3.1 The string branch. 

In the fiat case 7 = 0, equations admit the celebrated Nielsen-Olesen 

string as solution : 

N=l , L = X , P = Pno{x) , f = fno{x) (19) 

where P„o and /„o are determined numerically. When the parameter 7 is increased 
from zero, the solution (|T^ gets progressively deformed by gravity. In particular 
the functions A^, L become more complicated. Asymptotically they obey the 
following behaviour 

N{x ~^oo)=a (20) 

L{x ^ 00) = bx + c , b> (21) 

here a, b, c are constants depending on a, 7. According to (0), the Tolman mass 
vanishes. 

For fixed a, the set of solutions obtained by varying 7 (or vice versa) assemble 
into a branch of solutions called, after , the string branch. They exist as global 
solution on X G [0, 00] up to a critical value 7cr(«) (or, equivalently if 7 is fixed 
up to aeri'y)) which is reached when the parameter b (defined in (|2l|) becomes 
zero (solutions corresponding to 7 > 7cr- and 6 < will be the object of the next 
section). This phenomenon is illustrated by Fig.l for a = 1.0 (solid lines) and 
a = 3.0 (dashed lines). 

The figure further indicates that the inertial mass Mm increases slightly with 
7 while the parameter b decreases linearly with 7. Moreover, the critical value 
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7a(«) increases with 7, e.g. 

7,,(1.0) ^ 1.66 (22) 

7,,(2.0) = 2.0 (23) 

7e.(3.0) ^ 2.2 (24) 

The case a = 2 possesses further algebraical properties because there exist 
underlying self-dual (first order) equations implying the full equations (H), @. 
In particular we have in this case 

iV(x) = l. , M,„(7) = l (25) 

a=l , &=1-| , 7cr(2) = 2 (26) 

3.2 The Melvin branch. 

As discovered by , any global solution on the string branch possesses a shadow 
solution. This second set of global solutions of Eqs (§)-(§) assemble into a branch 
named, after the Melvin branch. The profile of a solution on the Melvin 
branch is presented in Fig. 2 for typical values of the parameters a = 1.8,7 = 1- 
In contrast to (|20|) , (p|), the asymptotic behaviour of the Melvin-branch so- 



lutions is such that 



N{x ^ 00) ~ Ax^/^ (27) 
L{x -> 00) ~ fix"^/^ (28) 

As a consequence the Tolman mass is not zero any longer : 

GMto = I hm {LNN') = \a^B (29) 

The inertial and Tolman masses of the string and Melvin branches are plotted 
on Fig. 3 as functions of a and for 7 = 1 (in this case the critical value of a is 
Uer ~ 0.155). The figure clearly indicates that both (inertial and Tolman) masses 



are higher for the Melvin branch than for the string branch. In the hmit a — Ocr, 
the numerical analysis confirms that the inertial massess (and also the Tolman 
masses) of the string and Melvin branches tend to a common value (in the case 
of the Tolman mass, this value is zero). 

The way the different radial functions associated with the two solutions ap- 
proach each other in the limit a — is illustrated by Figs 4 and 5 where the 
profiles of the functions are presented respectively for a = 0.2 and a = 0.16. In 
particular, it is seen on the figures that the matter functions /, P of the string 
branch deviate rather slowly from their conterparts of the Melvin branch. 

The evolution of the masses Mj„, Mfo for 7 ^ on the Melvin branch is sum- 
marized by Fig. 6. The radial integral defining these quantities clearly diverges 
when 7 approaches zero but the Einstein-Maxwell equations can be recovered 
from Eqs. (^-(§) by rescaling the radial variable x and the function L{x) ac- 
cording to 

x=^y , L = .I^L (30) 
V a 

and setting 7 = afterwards. Equation (P) then decouples and the remaining 
equations are the Einstein-Maxwell equations in the axially symmetric ansatz. 
The profile of the Melvin solution is illustrated on Fig. 7. It has 

GM.,n = \ (31) 

independently of a. 

This situation is somehow reminiscent to the case of the gravitating sphaleron 
solution of the Einstein- Weinberg-Salam equation (EWS) 0, [§]. The role played 
here by the flat Nielsen Olesen is played in EWS equation by the flat klinkhamer- 
Manton Sphaleron while the role of the Melvin solution on the top of the upper 
branch is played by the first solution of the Bartnik-McKinnon series of solutions 
0. However, unlike in the present theory, the transition from the sphaleron to 
the Bartnick-Mc Kinnon solutions in the EWS model is smooth Q]. 
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3.3 Solutions of vorticity > 1. 

Coming back to the flat Nielsen-Olesen equations, we know that solutions of 
vorticity n can be constructed while imposing the following boundary conditions 
on P{x) : 

P(0) = n , P(oo) = (32) 

Denoting by M{n) the classical energy of the Nielsen-Olesen solution of vor- 
ticity n, one of the distinguished feature of the self dual case (a = 2) is the mass 
relation 

M(n) = nM(l) (33) 

It is therefore a natural problem to check if the relation above still holds in 
presence of gravity. We did so by chosing 7 = 1 and by studying the inertial 
mass (which directly generalize the classical energy of the flat case) for the values 
of a close to a = 2. These results are presented on Fig. 8. The quantities Mj„(l) 
and |Mj„(2) cross at the self dual point a = 2 and their behaviour demonstrates 
that the Rebbi- Jacobs phenomenon also holds when gravity is added. Namely 
for a > 2 (resp. a < 2) the binding energy of the n = 2 string solution is negative 
(resp. positive). 

4 Closed solutions. 

As said above, the critical value 7cr(«) (or equivalently acri'y)) is determined when 
the slope parameter b of the function L si zero. For 7 > 7u(a) (or a < Quil)) 
the string and Melvin branches of solution continue to exist but they are not 
any longer available as regular solution for x G [0, 00]. The fact that one of the 
functions N 01 L has a zero at, say x = xq, leads to a singular point of Eq. (P) 
or d^). Correspondingly the solution can only be constructed for x G [0,Xo]. 

The inertial masses of the solutions continuing the Melvin and string branches 
for a < acr are presented in Fig. 3. Clearly, the closed solution having the lowest 
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mass is the continuation (for a < acr) of the Melvin branch. In it is called 
the Kasner branch. An example of such a solution is presented on Fig. 10 for 
7 = l,a = 0.14. It terminates at xq — 26. The evolution of xq, the maximal 
value of n, in function of a is presented on Fig. 11. 

5 Conclusion and outlook 

We have seen that, for a fixed a, the mass of a string grows with 7 while the 
Melvin-like configuration display the opposite behaviour. This results in the 
existence of a critical value ao such that the new solution is heavier (resp. lighter) 
than a string for a > ao (resp. a < ao). The numerical calculations show that 
tto = 0.5. Another particular value of interest is a = 2, for which there exists 
first order equations implying the classical equations (this is self duality). For 
a > 2 the Melvin-like configurations have a negative binding energy and, likely, 
are stable versus decaying into several solutions of unit vorticity. This is also the 
case with cosmic string solutions. So the kind of solutions considered here obey 
the Jacobs-Rebbi phenomenon on both sides of the self dual point a = 2. We 
also studied how the radius of the closed solutions grows with a. 

The Melvin-like solution raises many questions. A first one is a full under- 
standing of its geometry . A second one is the elaboration of a plausible mech- 
anism for the generation of such a configuration. For example, cosmic strings can 
be formed during a phase transition. This scenario seems conceptually difficult 
to implement here since one would have to perform quantum field theory at finite 
temperature and on a curved background. 

Cosmic string can also be formed by the collision of two monopoles. If the 
Melvin-type configuration has to be generated in this way, then the problem is 
a technical one since one has to solve the field equations for a fully time and 
position dependant ansatz. If any of these two mechanisms takes place with a 
reasonnable efficiency, then the new solution is likely to play a significant role in 
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cosmology, like its string cousins whose possible implications ranges from density 
generation to baryogenesis 

The kind of calculations reported here could be extended to other members 



of the Abelian models hierarchy ||10|| for which self dual equations are available as 
well. Lagrangians constructed by superposing a few members of the hierarchy 



could also be considered [|Tl|]; for such models, however, (and at least for flat 
space) only quasi self dual equations exist. 
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Figures Captions 

Figure 1 The inertial mass Mj„ and the quantities a, b defined in (^), (|2T|) as 
functions of 7 and for a = 1, (sohd) and a = 3.0 (dashed). 

Figure 2 The profiles of the solution on the Melvin branch for a = 1.8, 7 = 1.0. 

Figure 3 The intertial mass Mj„, the Tolman mass Hto and the quantities a, b 
of Eqs. (|20|) , (^TD as functions of a (for 7 = 1) for the string (solid) and 
Melvin (dashed) solutions. 

Figure 4 The string (solid) and Melvin (dashed) solutions for 7 = 1.0, a = 0.2. 

Figure 5 The string (solid) and Melvin (dashed) solutions for 7 = 1.0, a = 
0.16. 

Figure 6 The inertial mass Mj„ and Tolman mass Mto for the Melvin branches 
as functions of 7 for a = 1. The dashed lines represent the ratio M/7. 

Figure 7 The approach of the Melvin solution (7 = 0) by the solutions of the 
Melvin branch for decreasing values of 7. 

Figure 8 The quantity Mj„(n)/n for n = 1 (solid) and n = 2 (dashed) in 
function of a for 7 = 1. 

Figure 9 The profiles of gthe solution on the Kasner branch for 7 = 1, a = 0.14. 
The function crosses zero at x ~ 26. 

Figure 10 The evolution of the quantity Xmax on the Kasner branch in function 
of a and for 7 = 1. 
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